Flow-induced Agitations Create a Granular Fluid 
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We fluidize a granular medium through localized stirring and probe the mechanical response of 
quiescent regions far away from the main flow. In these regions the material behaves like a liquid: 
high-density probes sink, low-density probes float at the depth given by Archimedes' law, and drag 
forces on moving probes scale linearly with the velocity. The fluid-like character of the material is set 
by agitations generated in the stirred region, suggesting a non-local rheology: the relation between 
applied stress and observed strain rate in one location depends on the strain rate in another location. 

PACS numbers: 47.57.Gc, 83.80.Fg 



What governs the flow of granular media? At the grain 
level, interactions are mediated by collisions and contacts 
[lj|. While rapid flows where collisions dominate can be 
described by advanced kinetic theories Q, and the un- 
derstanding of flows where both contacts and collisions 
are important has recently advanced tremendously 0, 0| , 
it remains difficult to describe slow flows where enduring 
contacts dominate the interactions. 

Aspects of such slow grain flows can be captured by a 
frictional rheology in which the friction laws acting at the 
grain scale are translated to effective friction laws for the 
stresses acting at a coarse-grained level (H], [&, ?| ■ In such a 
Mohr-Coulomb picture, granular media remain jammed 
when the ratio of shear r to normal stresses P is below a 
critical value given by an effective friction coefficient, /x, 
while slowly flowing grains correspond to stresses close 
to the yielding criterion: t/P « fj,. 

This framework is, however, not complete. The com- 
bination of rate independence and a sharp yielding cri- 
terium leads to a description which predicts the local- 
ization of flows in shear bands of vanishing width and a 
corresponding sharp separation between stationary zones 
and flowing zones [Bj]. However, in experiments shear 
bands arc found to be of finite width and the boundary 
between flowing and stationary zones is not sharp, with 
creep flow occurring even far away from the main shear 
band @, @, EH EI El El • The first key question is there- 
fore: what is the nature of the nearly-stationary zones far 
away from the main flow? A second key question is mo- 
tivated by the observation that, for slow flows, the flow 
rate is independent of the stresses. But if the flow rate is 
not determined by the stresses what then is the physical 
mechanism that sets the flow rate of slow granular flows? 

Here we address these questions by locally stirring a 
container of glass beads at a rate f2 while probing the 
mechanical response of the essentially quiescent regions 
near the surface, away from the shear band (Fig. la- 
b). We find that a heavy object — such as a steel ball 
- sinks slowly into the sand with a rate proportional 
to £1 (Fig. la). Moreover, low density objects sink (or 
rise) to the depth predicted by Archimedes' law (Fig. lc). 




FIG. I: (Color online) (a) Snapshots of a stainless steel ball of 
diameter 25 mm and mass 64 g, sinking into a stationary gran- 
ular fluid. The granular fluid is generated in a split-bottom 
shear cell (inner radius R c = 81 mm, disk radius R s — 60 
mm) filled to a height of H = 60 mm with millimetric glass 
beads and driven at a rate Q — 0.1 rps. (b) Experimental 
setup, with the dome-like shear zone indicated in dark grey, 
(c) Examples of probe position, z, as function of time, t, for 
a filling height of H = 60 mm. A high-density object sinks in 
the sand (orange M = 48 g, D - 18 mm, !J = 3x 10~ 3 rps). 
A low-density object (M = 40 g, D = 40 mm, fi = 0.1 rps) 
sinks until reaching an equilibrium depth (red) — and rises 
up to this same depth if initially deeply submerged (blue). 



Therefore, granular materials do not exhibit a yield stress 
in the presence of flow: flow fluidizes granular media. By 
observing the motion of probes immersed in the sand, we 
find that the drag forces acting on the probes are linear 
in velocity. This suggest that the material is viscous, 
although we find that the viscosity depends strongly on 
filling height, location and mass of the probes — the 
rheology of the material is highly nonlinear. Moreover, 
the relation between applied stress and observed strain 
rate in one location depends on the strain rate in another 
15j. Our findings highlight novel behavior 
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which we believe to be crucial for the dev elop ment of 
better models of slow granular flows [l5, 1(| 17 1. 
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FIG. 2: (Color online) (a) Evolution of the probe position (z) 
for fi = 0.1 rps, H — 50 mm, D — 40 mm and five different 
masses ranging from 40 to 81 g. (b) The equilibrium depths, 
z eq , as a function of probe mass for H = 50,60 mm and 
D = 40 mm. The red diamonds and blue squares correspond 
to the sinking and rising cases as shown in panel (a). The 
green triangles correspond to lighter probes detached from 
the inductor measurement device and imaged with a CCD 
camera. The straight lines show Archimedes' law for a density 
of 1.92 kg/L and finite size correction factor k = 1.08 — the 
respective effective heights, H — SH, are 49 and 59 mm. 



Setup — Our experiment consists of a split-bottom 
shear cell [6, 11, ID, IH, filled with glass beads of mean 
diameter d = 1 mm to a height H (Fig. [T]) . The bottom 
disk (radius Rs = 60 mm) is driven by a micro-stepper 
motor at a rate, f2, which ranges from 1CP 4 to 1 rps. The 
relative humidity of the system is controlled at 8 ± 2 % at 
room temperature. Prior to beginning a measurement, 
the grains are stirred with a rod, the grain surface is 
leveled, and the bottom disk is spun at a rate of 0.5 
rps for 20 s. However, we have found no evidence for 
systematic dependence of the long-time behavior of the 
probe on the preparation history. Grooves and dimples 
are machined into the surface of the container and the 
disk to create a rough boundary and the resulting grain 
flows have been studied extensively @, [ll], [IH, [l6| ■ 

In order to investigate the liquid-like properties of the 
system, smooth aluminum cylinders of mass M and di- 
ameter D arc immersed in the grains. These probes are 
attached to a shaft that passes through an air-bearing, 
which fixes the horizontal position of the probe while al- 
lowing the probe to rotate and to move in the vertical 
direction. The probe shaft next passes into a LVDT po- 
sition sensor (DC Fastar 2M, accuracy of 2 fxm) which 
we employ to measure the position of the bottom surface 
of the probe, z. 

Phenomenology For the large filling heights 
(H/R s > 0.8) examined here, the shear zone that em- 
anates from the edges of the disk is submerged in the 
bulk of the material (Fig. lb), and the residual flow ob- 
served at the free surface is several orders of magnitude 
smaller than the driving rate 11 , I2| ■ A probe placed 



and M — 100 gr, the initial depth is less than 5 mm). 
However, when the disk begins rotating the probe will 
immediately start to sink into the grains. Heavy probes 
will continue to sink until they reach the bottom of the 
container, while light objects reach a well-defined float- 
ing depth (Fig. [TJ;). By observing the motion of such 
probes we will address the following questions: What sets 
the floating depth of light objects? Does the material ex- 
hibit a yield stress? What governs the drag forces on the 
probes? 

Archimedes' law — As shown in Fig. [5^, we have de- 
termined the floating depths (z eq ) for a series of probes 
of varying mass. For each mass, three sinking and three 
rising experiments are conducted — in all cases the equi- 
librium depths, z eq , are consistent within 0.6 mm. We 
conclude that there is no appreciable yield stress in the 
material: if the system did exhibit a yield stress, the sink- 
ing and rising probes would reach different equilibrium 
positions. As shown in Fig. [2]b, the equilibrium depth 
varies linearly with M. Both this linear dependence and 
the convergence of the rising and sinking probes to the 
same floating depth are suggestive of Archimedes' law, 
which in this case reads: 



z eq = H(z eq )-SH - 



4M 
n(D) 2 p 



= H -SH- 



AM 
irD 2 np 



(1) 



at the surface of the beads will get stuck at a small depth 
if the disk is not rotating (for example, for D = 30 mm 



Here, both D and k are effective parameters that incorpo- 
rate finite size effects. First, the probe displaces grains as 
it sinks, so the actual height of the grains, H(z), depends 
on the probe position. Second, the finite size of the grains 
suggests that the effective size of the probe, D, may be 
somewhat larger than the real value: D < D < D + d. 
Both effects can be incorporated in the finite size param- 
eter k — the first effect is of order l + (D/(2R c )) 2 w 1.06, 
the second effect is of order (D/D) 2 which ranges from 
1 to 1.05 — hence we estimate k to be between 1.06 and 
1.11. Finally, we allow for a small correction, SH, which 
is on the order of a grain size and takes into account 
the lower packing density of grains near the bottom of 
the probe. Direct measurement of the mass and volume 
of well-compacted grains yields a density of 1.92 ± 0.05 
kg/L. 

As shown in Fig. ^jp, the equilibrium depth is well de- 
scribed by equation |T]) for n = 1.08 and SH = d, where d 
is the diameter of the beads. The finite size correction, re, 
lies within the expected range, and the height correction 
SH is also small — we thus conclude that Archimedes' 
law describes the floating depths of our probes accurately, 
provided that finite size corrections are properly taken 
into account. Note that Archimedes' law has also been 
observed in a granular system in which the lateral bound- 
aries are vibrated [i~8| — but this driving appears far 
more vigorous than in our system. 

Viscous Force — As shown in Fig. [3J the probe position 
approaches the equilibrium depth exponentially for small 
\z — z eq \ The data suggests that the characteristic 
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FIG. 3: (Color online) Exponential approach to equilibrium 
depth as demonstrated by sinking probes (red crosses) and 
rising probes (blue circles) for M = 40 g, D = 40 mm, fi = 0.1 
rps and H = 50 mm (a) and H = 60 mm (b). 



times for rising probes are somewhat smaller than that 
for sinking probes, although the apparent slope of the 
curves is very sensitive to the precise choice of z eq . In 
the remainder of this paper we focus primarily on sinking 
probes where we can vary the probe mass over a wide 
range. 

Since the difference between the gravitational and 
buoyant forces is linear in z — z eqi the exponential re- 
laxation implies that the drag force (F^) on a probe is 
proportional to its velocity. Ignoring geometrical form 
factors for drag on a cylinder, we define the effective vis- 
cosity, t], of the granular fluid using Fd = —rjDdz/dt. 
Combining this drag force with the buoyant force, F bl 
and the gravitational force suggested by Archimedes' law 
yields the following equation of motion for the probes: 



dz 

1] D— = -Mg + F b 
dt 



(2) 



In contrast, drag forces on intruders in quiescent or os- 
cillated granular media reported before 20, 21 1 exhibit a 
finite threshold. 

In the remainder of this paper we will address the fol- 
lowing two questions, (i) How does the effective viscosity 
depend on the parameters H and Ml (ii) Is the viscosity 
set by the local residual flow near the probe? 

Viscosity — From a wide range of experimental data 
we conclude that the viscosity is inversely proportional to 
Vt. In Fig. 2^, we plot the immersion speed, v, at a fixed 
depth (z = H—13 mm) as a function of £1 for a number of 
probes. Clearly, v oc Q over a wide range of disk rotation 
speeds. The characteristic time scale of the exponential 
relaxation of floating probes is also proportional to Q (not 
shown). Hence, the relevant timescale for probe motion 
is set by f2, and internal timescales (such as given by 
vibrations) appear irrelevant. 

In a true liquid, the viscosity is independent of the 
mass of an object sinking in the system. However, the 
viscosities of our granular liquids exhibit a surprising de- 
pendence on the probe mass (Fig. HJa). Moreover, we 



FIG. 4: (Color online) (a) Probe velocity 13 mm below the 
surface as a function of disk rotation speed Q. Red stars: 
H = 60 mm, M = 77 g, D = 18 mm, Blue diamonds: H = 50 
mm, M = 46 g, D — 18 mm, Purple squares: H = 60 mm, 
M = 48 g, D = 18 mm, Green crosses: H = QG mm, M = 86 g, 
D — 18 mm. Orange triangles: H = 60 mm, M = 46 g, D = 18 
mm with a probe immersed at a radius of 15 mm away from 
the center. The line has slope 1. (b) Rescaled viscosity, £ln, 
at z = 50 mm as a function of probe mass for D = 40 mm, 
H = 60 mm. The line is a guide to the eye corresponding to a 
power law with exponent -1.5. 



found that the precise form of T)(M) depends strongly 
on H and weakly on the measurement depth, z. The 
detailed rheology is complicated and we leave a detailed 
study to further work. 

In Fig. [3Ji we plot the measured viscosity for a single 
probe (M = 59 g, D = 30 mm) as a function of z for three 
values of H (50, 60 and 70 mm). These viscosities are 
obtained by measuring z(t) and then invoking Eq. @ to 
calculate the effective viscosity. From Fig. [5jt it is appar- 
ent that the viscosity depends on the filling height, H . 
Moreover, the viscosity for a given filling height changes 
very little with the probe depth, which is consistent with 
our observation that the probe position of light probes 
evolves exponentially over a substantial range. Clearly, 
the strong dependence on H and the weak dependence 
on z rules out a picture in which the viscosity simply 
depends on the distance to the flowing zone. 

Nonlocal Flow Rule — Even though the flow and strain 
rates near the free surface are very small, they are not 
zero [f| [ll], [l3] ■ It is therefore instructive to ask if the 
probe motion is determined by the local residual flow 
near the probe. To answer this question, we will compare 
the H and z dependence of the viscosity, shown in Fig.[5^, 
to the dominant local strain rate. For dome-like shear 
bands, the strain in the region above the center of the disk 
is torsional 12| , and so the strain rate varies as {d z uj\ r —o). 

We have measured the precession rate u> p (z) := 
u)(z)\ r= o/tt by inserting tiny vane-like probes in the cen- 
ter of the grain flow and observing the rotation of the 
probe with a rheometer. The resulting precession rates, 
u) P (z), are shown in Fig. [5b for H = 50,60 and 70 
mm, and are in good qualitative agreement with ear- 
lier MRI measurements and simulations Our larger 
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FIG. 5: (Color online) (a) Viscosity as a function of probe 
position for a probe with M — 59 g, D — 30 mm and H = 50, 60 
and 70 mm (lower orange curve, middle red curve and upper 
blue curve respectively), (b) The variation of the strain rate 
d z uj p with z for H = 50, 60 and 70 mm. Inset: ui p as a function 
of the height above the spinning disk. Curves are of the form 
lo p (z) - lu p (z = H) ~ exp(-(z/^) 1 ' 5 ), with £ = 9.5, 11.5 and 
13 mm respectively. 



measurement range allows us to establish that u) p (z) ap- 
proaches the surface precession rate, lo p (z = H), faster 
than an exponential but slower than a Gaussian, and 
our results are fitted well by an expression of the form 
u)p(z) — Lu p (z = H) = id p (z = 0) cxp( — ), where 

lu p (z = 0) captures slip near the bottom disc, and the 
characteristic length scale £ is of order 10 mm. 

By either differentiating this expression, or by numeri- 
cally differentiating the measured data, we can determine 
d z oj p (z): the result is shown in Fig. [5)3. Clearly the local 
strain rate and r){z) correlate very poorly — while the 
former changes over four decades for H = 70 mm, the 
latter changes over less than half a decade. 

Conclusion and Outlook — By applying localized shear 
to a container of glass beads we have created a granular 
fluid with unusual properties. The first surprise is that 
our granular fluid does not exhibit the glassy behavior 
typical of agitated granular media [H, |23, 0, HH > and 
that the characteristic time scale is simply set by the in- 
verse driving rate I/O. We speculate that the application 
of steady shear does not allow for caging and other typ- 
ical glass-like effects, but instead leads to fluidization of 
the material. 

The second surprise is that the drag forces on the in- 
truder are simply proportional to the velocity, which is 
in contrast to earlier work [2^, 21 1. The third surprise 



is that our granular fluid does not exhibit a finite yield 
stress. The disappearance of the yield stress in the pres- 
ence of flow suggests that the stress and strain tensor are 
co-linear, consistent with recent simulations (ljj . 

The final surprise is that the local rheology of the ma- 
terial is set by the flow in the stirred region, and not by 
the local, residual (creep) flow. We suggest that random 



grain motion in the flowing zone leads to grain agitations 
even far away from the flow, and that these agitations 
generate the liquid-like behavior. Hence, the rheology 
is nonlocal: the relation between applied stress and ob- 
served strain rate in one location depends on the strain 
rate in another location. Whether this is similar to obser- 
vations of nonlocal rheology recently observed in emul- 
sions [3] is at present an open question. 
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